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We show that various factors can lead to an isotopic dependence of the penetration depth 8. Non- 
adiabaticity (Jahn- Teller crossing) leads to the isotope effect of the charge carrier concentration n 
and, consequently, of 8 in doped superconductors such as the cuprates. A general equation relating 
the isotope coefficients of T c and of <5 is presented for London superconductors. We further show 
that the presence of magnetic impurities or a proximity contact also lead to an isotopic dependence 
of 8; the isotope coefficient turns out to be temperature dependent, (3(T), in these cases. The 
existence of the isotope effect for the penetration depth is predicted for conventional as well as for 
high-temperature superconductors. Various experiments are proposed and/or discussed. 



I. INTRODUCTION 



When speaking of the isotope effect (IE) in supercon- 
ductors, one generally considers the influence of isotopic 
substitution on the superconducting critical temperature 
T c . In systems where the pairing mechanism is at least 
partly mediated by the electron-phonon interaction T c 
depends on phonon energies. Thus, replacing one ele- 
ment by its isotope will affect T c via the. change in phonon 
frequency. In a series of recent paperaJTJ we have shown 
that besides this most simple case, there are a number of 
other factors, not related to the pairing mechanism, that 
can strongly affect the isotope coefficient (IC) of T c . The 
influence of the proximity effect, the presence of magnetic 
impurities and non-adiabaticity were specifically studied. 
Several experiments were proposed for conventional su- 
perconductors and the oxygen IE in high-T c materials 
was discussed. In the latter case it was experimentally 
shown (see e.g. Refs. ||-|8| and references therein), that 
the value of the isotope coefficient of T c ranges from al- 
most zero to values above 0.5, depending on the doping 
process. The theory developed in Ref. g allowed us to 
present a systematic description of all experimental data 
collected to this day on the oxygen isotope effect of T c in 
YBa2Cu307_5 related systems. We stress here the fact, 
however, that the theory is not restricted to high-T c ma- 
terials, but applies to conventional superconductors as 
well. 

The influence orutJac IE of factors not related to the 
pairing mechanismETB raises the question, whether be- 
sides T c , there are other superconducting properties that 
also display these new effects. In the present paper we 
focus on one of the fundamental parameters, the penetra- 
tion depth, and demonstrate that, indeed, it can display 
a non-trivial dependence on isotopic substitution. The 
treatment is based on our previous analysis in Ref. |^ that 
focused on the unconventional isotopic dependence of T c . 

After isotope substitution in a superconductor [M — » 



M* = M + AM, M is the ionic mass), the value of the 
critical temperature is shifted: T* Q — T c0 + AT c0 . As is 
known, the isotope coefficient ao is defined by the relation 
T c0 ~ M~ a ° and is equal to a = -(M/AM)(AT c0 /T c0 ). 
Suppose now that the superconductor is perturbed by 
some external factor as e.g. the presence of magnetic 
impurities or a normal film (proximity effect). Super- 
conducting properties are affected by these factors. For 
example, T c q — > T c and, as a consequence, the depen- 
dency of the new critical temperature on the ionic mass 
(T c ~ M~ a ) will also be changed (a ^ a , see Refs. |,|). 
This not only affects the isotope effect of T c but also, e.g., 
the IE of the penetration depth. 

In the following, we are concerned with the isotopic 
dependence of the penetration depth. By analogy with 
T c , we define the new isotope coefficient (3 by the relation 
S ~ M" 13 . Therefore, (3 is determined by 



M AS 



(i) 



where AS is the shift of the penetration depth induced 
by isotopic substitution. 

The structure of the paper is as follows. Section |n] is 
concerned with the non-adiabatic IE, whereas Sees, fill 



and [V address the impact of magnetic scattering and 
the proximity effect on the IC, respectively. An interest- 
ing property of the two latter factors affecting the IE is 
that the isotope coefficient appears to be temperature de- 
pendent. We not only demonstrate here that such effects 
exist but that they are also sizeable. Moreover, we ar- 
gue that the measure of the IC can be used to probe the 
presence of non-adiabaticity, magnetic impurities or the 
proximity effect in superconductors. The theory leads 
us to propose several experiments for conventional and 
high-temperature superconductors. We also discuss ha 1 
recent experiments done on high-T c superconductors!! 
can be described by our theory. 



II. ISOTOPE EFFECT OF S AND 
NON-ADIABATICITY 

The theory for the non-adiabatic IE of T c has recently 
been proposed^ and applied to describe the oxygen IC in 
high-temperature .superconductorsQu as well as proper- 
ties of manganitesEJ. We show in this section that other 
quantities like the penetration depth 6 of a magnetic field 
can be affected by isotopic substitution through this non- 
adiabatic channel. 



A. Theory 

1. Isotope effect of n 

Let us describe the main concept of the non-adiabatic 
IE introduced in Rcf. [j] (see also Ref. |J). Consider a 
system composed of a conducting system and a charge 
reservoir (as for example the Cu02 plane and the buffer 
layer of a cuprate high-T c superconductor) and assume 
that there is a charge transfer between these subsystems. 
In general, the reservoir-conducting-layer charge-transfer 
process occurs through a group of atoms located in both 
subsystems and/or bridging the two. For example, the 
charge transfer between the Cu02 plane and the reservoir 
(CuO chain for YBCO, BiO plane for Bi-based oxides, 
etc..) occurs through the apical oxygen. The electronic 
degeneracy and the corresponding Jahn- Teller instabil- 
ity of the group of atoms leads to the crossing of elec- 
tronic terms. Because of the breakdown of the adiabatic 
approximation for these atoms, the electronic and ionic 
degrees of freedom can no longer be considered as decou- 
pled. The potential energy surface considered as a func- 
tion of the ionic coordinate is composed of two electronic 
terms with two close minima. For example, the motion of 
the apical oxygen which is a lighter ion of the Cu-O-Cu 
Jahn- Teller complex is characterized by two minima (a 
more detailed discussion is given in Refs. 0J| ; experimen- 
tal evidence of these minima can be found in Refs. |IT| . |l2| ). 
Note that each of these terms can be considered in the 
harmonic approximation. The "double-well" is thus not 
related to the anharmonicity of the lattice, but to the 
non-adiabatic behaviour of the group of ions involved in 
the charge-transfer process. In a qualitative picture, the 
charge transfer between the conducting and the reservoir 
subsystems will involve the tunneling of the ionic com- 
plex between the two minima and, as a consequence, the 
density of charge carriers in the conducting subsystem 
will depend on the ionic masses. Properties depending 
on the charge-carrier density n will thus be ionic mass 
dependent (n = n(M)) and exhibit an isotope effect. 

In order to consider the Jahn- Teller crossing (where 
the adiabatic approximation is not applicable) it is con- 
venient to use the so-called diabatic representation, (see, 
e.g., Ref. |l4|). In this representation one can shown that 
the energy level splitting has the form 
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with 



L = f dripl (r, R )# e </>2(r, R ) 

F12 = /"dR^(R)$2(R) =Fi 2 (M) . (3) 

H e is the electronic part of the Hamiltonian and 
*j(r,R) = ipj(r, R)<pj(R) (j = 1,2). ^(r,R) is the 
electronic wave function (with electronic coordinates r) 
depending parametrically on the nuclear coordinates R 
and $j(R) is the vibrational wave functions. The last 
equality in Eq. (||) was obtained because ipj (j = 1, 2) is 
a slow varying function of R and can thus be evaluated at 
R , the crossing of electronic terms. L does not depend 
on the ionic masses. On the other hand, the important 
Franck-Condon factor F±2 depends on the lattice wave 
functions 3?j and thus on ionic masses M. 

The reservoir-conducting-layer charge transfer is ac- 
companied by a change of electronic terms and can be 
visualized as a multi-step process. First, the charge car- 
rier can move from the reservoir to the group of ions 
(e.g. from the chains to the apical oxygen in YBCO). 
Then, the complex tunnels to the other electronic term 
("Ji — > ^2)- As a final step, the charge carrier can hop 
to the conducting layer. We emphasize that this charge 
transfer is a dynamical process and leads to a stationary 
state in the sense that the charges oscillate in time be- 
tween the reservoir and the conducting subsystems, with 
a frequency given by Eq. (Q). Since the ionic masses affect 
F12, and thus ii/12, different isotopic masses imply differ- 
ent values of the characteristic charge-transfer frequency. 
This in turn affects the charge-carrier concentration n. 

The isotope coefficient of T c can be written as a — 
a ph + a na , where a ph = (M/T c )(dT c /dQ){dn/dM) is 
the usual (BCS) phonon contribution (f2 is a character- 
istic phonon energy) and the non-adiabatic contribution 
is given by: 



n dT c 



(4) 



where the parameter 7 = —M/n(dn/dM) has a weak 
logarithmic dependence on M (see Ref. |l|) . Eq. (^) shows 
that the IC of T c depends on the doping of the conduct- 
ing layer and on the relation T c {n). This result was used 
in Ref. [2] to analyse the IE of high-temperature super- 
conductors. 



2. Non-adiabatic isotope effect of 8 

Let us now focus on the non-adiabatic IE of the pene- 
tration depth. As shown above (see also Refs. non- 
adiabaticity for doped materials such as cuprates, leads 
to the dependence of the charge-carrier concentration n 



2 



on the ionic mass M, that is n = n(M). In the weak- 
coupling London limit the penetration depth is given by 
the well-known relation 



S 2 = 



47m s e 2 47rn(/?(T/T c )e 2 



(5) 



where m is the effective mass. n s is the superconducting 
density of charge carriers, related to the normal density 
n through n s — nip(T/T c ). The function ip(T/T c ) is a 
universal function of T/T c . For example, for conventional 
superconductors, ip ~ 1 — (T/T c ) 4 near T c , whereas tp ~ 1 
near T = (in the absence of magnetic impurities; their 
influence is discussed at the end of this and in the fol- 
lowing sections). Note that the relation 5~ 2 ~ n s is also 
valid in the strong-coupling case (see, e.g., Refs. ||). 

We can now determine the isotope coefficient (3 of the 
penetration depth from the relation (0): 



(3=- 



M dS dn„ 



M dn. 



5 dn, dM 2n, dM 
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Because of the relation n s — nip(T), one has to distin- 
guish two contributions to (3. There is a usual (BCS) 
contribution, [3 p h, arising from the fact that tp(T/T c ) de- 
pends on ionic mass through the dependency of T c on the 
characteristic phonon frequency. Indeed, isotopic substi- 
tution leads to a shift in T c and thus in 5, whi ch m ight 
be noticeable near T c (see the discussion in Sec. [IB). In 
this paper, however, we focus on the non-trivial manifes- 
tation of isotopic substitution, arising from the isotope 
dependence of the charge-carrier concentration n. 

From Eq. (||) and the relation n s — nip(T) it follows 
that 



P — flna + f3ph 



where 
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1 ~ 2ridM 
M dip{T) 

2tp(T) dM 
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(8) 
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Note that n(M) is the normaZ-state charge-carrier 
concentration. Thus, unlike f3 p h (which depends on 
<j>(T/T c )), the non-adiabatic coefficient (3 na does not de- 
pend on parameters characterizing the superconducting 
state. This effect should be observed near T — OK, where 
fiph is negligibly small. 

Comparing Eq. (Q) and Eq. (j^) one infers that j3 na — 
—7/2 and thus establish a relation between the non- 
adiabatic isotope coefficients of T c and 5: 



n dT c 



(10) 



This result holds for London superconductors. The equa- 
tion contains only measurable quantities and can thus be 
verified experimentally. It is interesting to note that (3 na 



and a na have opposite signs when dT c /dn > (which 
corresponds to the underdoped region in high-T c materi- 
als). 

In later sections we will discuss the influence of mag- 
netic impurities on the isotope effect of T c and 6. We note 
here that in the presence of magnetic impurities the re- 
lation n s — rup(T) remains valid, but tp(T) depends now 
on the direct scattering amplitude T2 defined in Sec. Ill B 
(this results, e.g., in the inequality n s (T = 0) < n in the 
gapless regime). As a consequence, magnetic impurities 
affect (3 p h in Eq. (Q), but leave (3 na and thus Eq. jlfj| ) 
unchanged. 



B. Non-adiabatic isotope effect in high-T c 
superconductors 

The value of the parameter 7 (see Eq. (|])) has been 
evaluated in Ref. || for Pr-substituted (YPrBCO) and 
oxygen-depleted (YBCO) YBa2Cu30y_5 and was found 
to be 7 = 0.16 and 7 = 0.28 respectively. With 
Eq. ( |To| ) one then obtains the IC of the penetration depth 
/3 na = —0.08 in the first case and (3 na = —0.14 for the 
latter. n 

Experiments-have been performed on YBCOu and 
La2- 2; Sr a ;Cu04EI (LSCO) which indicate that the Meiss- 
ner fraction (and thus the penetration depth) indeed dis- 
plays an isotope shift. Unfortunately, the measurements 
on YBCO have been done near T c where any supercon- 
ductor (conventional or bigh-T c ) displays an isotope shift 
of 5 through the dependency <S(T/T C ) - y/l -T 4 /T c 4 and 
corresponds thus to the usual BCS IE (this is the con- 
tribution (3 p h arising through (p(T) as discussed in the 
last section; see also Ref. |l5|). The BCS isotope coeffi- 
cient of S (near T c ) can easily reach values of the order 
observed in the experiments, even for a very small IC of 
T c (with a value a = 0.025 as observed in YBa2Cu307_s 

one obtains (3 0.6 for T/T c - 0.95). To avoid the 

BCS contribution to the IC that appears in all supercon- 
ductors because of the strong dependency S(T/T C ) near 
T c it is better to study experimentally the IE of the pen- 
etration depth near T = 0. This statement is valid if one 
is interested in determining the non-adiabatic contribu- 
tion to the IC (3 but is not general. As will be shown 
below, it is for example possible to extract the influence 
of magnetic impurities on the IE even near T c . 

The situation with the. measurements of Ref. |^ is dif- 
ferent. The experimentsEI done on LSCO were obtained 
near T = and reflect the unconventional dependence 
S(M) (the BCS contribution to the IC of 6 vanishes at 
these temperatures). The shift has been measured near 
x ~ 0.11 and x ~ 0.15. 

Near the concentration x ~ 0.11, T c experiences a 
small depression as a function of doping x&B Although 
several explanations have been proposed (as, e.g., the 
presence of electronic inhomogeneities, a change of the 
electronic density of states due to lattice distortions, im- 
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purity scattering, magnetic ordering, etc... 0) the pres- 
ence of the dip is not well understood yet. The choice of 
this concentration is thus unfortunate and inappropriate 
for the study of the IE. It was arguedB that there is no in- 
fluence of isotopic substitution on the charge carrier con- 
centration n near x ~ 0.11 and that n can therefore not 
depend on the ionic mass. To support this idea the influ- 
ence of isotopic substitution on the structural (tetragonal 
to orthorhombic) transition temperature T s was studied!! 
The assertion is based on the assumption that T s depends 
solely on n, that there is a monotonic relation between 
T s and n with non-zero slope and a one-to-one corre- 
spondance between p-,and x. However, the experiments 
performed on LSCOtZl show that the hole concentration 
n dependency of T c and T s as well as the relation n(x) 
are not well established, especially in the region of the 
dip. The structural phase transition temperature T s has 
also only beeij.-pieasured at three points in the vicinity 
of the dipEj'EJB. Since T c has an unexpected behaviour 
it would be of interest to study the detailed dependency 
T s (x) in this region. Furthermore, as known from other 
high-temperature superconductors (as, e.g., YBCO) it is 
unlikely that the relation between n and x is linear. We 
emphasize that our analysis of the IC from Eqs. (0) and 
( |l0| ) does not require any assumption on the relation n(x) 
since we rely solely on the experimental relation as ob- 
tained, e.g., from /iSR experiments (in this case one has 
T c (<7 ~ n s /m*); see Ref. yj). In short, since the structural 
and electronic properties are not well understood in the 
region around x ~ 0.11, it does not allow a conclusive 
statement on the IE of 5. n 

The only experimental resulto that can be discussed, 
and that does not contain the BCS contribution is the 
one done at x ~ 0.15 (for LSCO). Using Eq. (0) and 
the value ofr-,the experimentally observed relative shift 
A5/5 = 2%0 one obtains /3 na ~ 0.16 [and 7 = 0.32, 
cf. Eq. and (|ic|)1 which is in good agreement with the 
calculations presented above. Indeed, it is of the same or- 
der as (3 na for oxygen-depleted YBCO. This is consistent 
with the fact that in both these materials the reservoir- 
Cu02-plane charge transfer involves the same group of 
ions including the apical oxygen. The larger value for 
LSCO can be traced back to the fact that in YBCO one 
has taken into consideration the presence of magnetic im- 
purities (which also enhances the IC), whereas in LSCO 
there is no indication of rfche presence of significant quan- 
tities of such impurities fl c 

One notes further that the IC of 5 is much larger in 
YBCO and LSCO than in YPrBCO. The reason for this 
discrepancy can be twofold. On the one hand, the charge- 
transfer channel is different in YPrBCO from the two 
other materialso (e.g. it does not involve the apical oxy- 
gen). On the other hand, the concentration of magnetic 
impurities is higher in YPrBCO than in YBCO (and 
LSCO). This was taken into account in the calculation 
of a na and influences the value of (3 na in YBCO related 
systems as well. A larger part of the IC is thus due to 
the magnetic impurities in YPrBCO (see Ref. [2] for a 



detailed study of the interplay between non-adiabaticity 
and magnetic impurities). 

One should add that the experiment determining|-jthe 
IC through magnetic susceptibility measurementua is 
rather inaccurate. More reliable data could be ob- 
tained by measuring the penetration depth shift using mi- 
crowave measurements or the Josephson effect. It would 
be interesting to perform experiments on the systems 
discussed above that display the non-adiabatic charge- 
transfer channel. 

As stressed earlier and according to the analysis of 
Ref. U one has, however, to include the contribution of 
magnetic impurities to describe correctly the IE of T c 
in some of the high-T c superconductors. The magnetic 
impurities also directly affect the IC of 5 and should be 
treated as well. This will be done in the next section. 



III. ISOTOPE EFFECT OF 5 AND MAGNETIC 
IMPURITIES 

As described in Refs. ||J^ magnetic impurities modify 
the IC of T c . We show here that it can also induce an 
IE of the penetration depth. Consider a superconductor 
(conventional or highT c Ldoped with magnetic impurities. 
Abrikosov and Gor'kovEj have shown that these impuri- 
ties act as pair breakers. In the dirty limit one can calcu- 
late the penetration depth analytically for temperatures 
near T c and at T = 0. 



A. Magnetic impurity contribution to (3 near T c 

The penetration depth 5 in the presence of magnetic 
impurities and near T c has been calculated by Skalski et 



al 
by 



Taking terms up to the order A the result is given 



(ii) 



where a = Aajy/c (<7n is the normal state conductivity), 
C(z. <l) = £„>o l/(«-+<?) 2 and 7,, = r s /27rT c (r s = f s n M 
is the spin-flip scattering amplitude; iim is the concentra- 
tion of magnetic impurities and F s is a constant). Near 
T c , the order parameter can be written asE3: 



A 2 = 2F 2 (1 -r) 



1 - C2 + (1 - t) [1 - & + C" 3 ] 
C3-C4 



V* (1= 



with C, = 7f _1 C(2, Is + 1/2), z = 1, 2, ... and r = T/T c . 
Using Eqs. (0) and @ one can calculate the isotope co- 
efficient of the penetration depth in the presence of mag- 
netic impurities (/3 m ). In order to single out the impact 
of magnetic impurities on the value of the IC, we calcu- 
late the quantity (3 m = -(M/AM)(A5/5) [cf. Eq. 
where 5 = 5(T, T S )/5(T, 0). As a result, the IC is written 
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as a difference (3 m = m — {3o where (3 m (0o) is the iso- 
tope coefficient of S(T,T S ) (8(T,0)). Using the relation 
(d/dq)C(z } q) = —z((z + 1, q) one can express (3 m in the 
form: 



Pm = (Ri — Ra)ot r , 



where 
Ri = 
Ro = 



-Ji 



N2 
Ni 



£2 



(2 



1 



~\— /o = ~ [1-^(7. + 1/2)7.] /o 

2 Ctm 2 



(13) 

(14) 
(15) 



and f = (3 - r 2 )/(!_- r)(3 - r). The functions iVi, L>i 
are defined in Eq. ( Jl2| ) and 

Y 2 = 3(1 - r) 2 [C 2 - 2C 3 + U] + t{2 - t) - & 
£> 2 = 2 (3C 4 - Ca - 2C 5 ) 

Furthermore, Eq. ( |l5| ) has been written using the rela- 
tion for the isotope xaefHcient of T c in the presence of 
magnetic impurities oEj 



Of) 



1 - ^(7. + l/2) 7s 



(16) 



where ao is the isotope coefficient of T c0 , the critical tem- 
perature in the absence of magnetic impurities. -0' is the 
derivative of the psi function. 

Eq. ( |l3| ) shows that the isotope coefficient of S is pro- 
portional to the isotope coefficient of T c . This relation is 
valid near T c (where A is small) and for impurity concen- 
trations such that the condition AT C /T C <gc 1 is satisfied. 

An important feature of the result (|l3|) is that the IC 
of 6 is a universal function (3(T/T C ). Indeed, Eq. @ 
can be calculated once T/T c , T c0 , ao and 7. are given. 
The three first quantities are obtained from experiment. 
The last parameter 7. also determines the depression of 
T c induced by magnetic impurities antLcan be calculated 
from the Abrikosov-Gor'kov equationLJ: 



In 
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V> o + 7e 







(17) 



From Eqs. @, @ and @ we thus conclude that the 
knowledge of the measurable quantities ao, T c0 and T/T c 
completely determines the IC f3 m . 

If besides the magnetic impurity channel there is also 
a contribution from the non-adiabatic channel as pre- 
sented in the last section, the total IC takes the form 
Ptot = Pph + (3 m + Pna- In addition, according to Ref. |, 
one has to replace ao in Eq. ( |l6| ) by a na + a p h (a na is 
given m Eq. (§ and a p h is the usual phonon contribu- 
tion to the IC of T c ). The non-adiabatic and magnetic- 
impurity contributions to the IC are thus non-additive 
near T c (because of the presence of a na in (3 m and the 
presence of T2 in j3 p h] as shown in Sec. |II A 2 
fected by magnetic impurities) 



_ (i na is unaf- 
One notes further, that 



if the IC is induced by magnetic impurities only, then its 
value is always negative, whereas the situation might be 
different in the presence of a non-adiabatic contribution. 

The factor /3q in Eq. (|l^) which is the IC of S in the ab- 
sence of magnetic impurities contains a trivial tempera- 
ture dependence due to the factor <p(T/T c ) - 1 - (T/T c ) 4 
and the isotopic dependence of T c . However, the presence 
of magnetic impurities leads to an additional, novel tem- 
perature dependence of /3 m , beyond the BCS one. This 
temperature dependence is shown in Fig. 1 for fixed val- 
ues of T c (i.e. of the magnetic impurity concentration). 
The solid and dashed lines describe the case where only 
magnetic impurities are added to the system (no non- 
adiabatic contribution). The situation corresponds to 
Zn-doped YBCO (YBCZnO). It has been shown with use 
of NMR technique that Zn-substitution (for Cu) occurs 
in the Cu02 plane and leads to the appearence of local- 
ized magnetic moments (see Ref. ^l] and below). The two 
other curves (dotted and dash-dotted lines) are obtained 
when, in addition, a non-adiabatic channel is present. 
They are discussed below. Fig. 1 shows that the absolute 
value of the IC increases as one approaches T c . This tem- 
perature dependency is significant since in certain cases 
the value of (3 triples when r = 0.75 increases to r = 0.95. 
Note also, that the solid and dashed lines describe con- 
ventional superconductors as well, since only ao has to 
be modified (in YBCZnO ao = a p h = 0.025 which is 
much smaller than typical values for conventional super- 
conductors). 

Fig. 2 shows the dependency of the IC j3 m on T c (i.e. on 
the concentration tim) for fixed values of the tempera- 
ture. It appears that an increase of magnetic impurity 
concentration leads to an increase of \f3 m \- Moreover, the 
change of (3 m (T c ) is stronger as T approaches T c . It turns 
out that the qualitative behaviour of (3 m (T c ) is similar to 
ot m {T c ) [Eq. ( [lq ) with on = ot p h] but with opposite sign 
(see also Fig. 1 of Ref. ||) . As for Fig. 1 , the parameter 
ao corresponds to the value for YBCZnO, but the quali- 
tative picture holds for conventional superconductors as 
well. 

The calculation of the IC near T c presented above sug- 
gests an experiment to determine the contribution of 
magnetic impurities to the IC. For conventional super- 
conductors where the magnetic moments are introduced 
into the superconductor as impurities one can substract 
the BCS isotope coefficient discussed at the end of last 
section by measuring its temperature dependency. For 
high-temperature superconductors the same procedure 
can in principle be applied. Nevertheless, one has to be 
more carefull since magnetic impurities are intrinsically 
present in some of these materials as shown in Ref. |2^. 
Indeed, there are evidences that already optimally doped 
YBCO contains magnetic impurities. Upon depletion of 
oxygen one further adds magnetic impurities. One can 
thus not avoid the presence of these localized magnetic 
moments. As a consequence, it is more difficult to extract 
the BCS component for these materials. The procedure 
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in high-T c superconductors would thus consist in measur- 
ing the (temperature dependent) isotope coefficient of 5 
at optimal doping (which contains the smallest concen- 
tration of intrinsic magnetic moments) and substract this 
component from the measurements at a given doping. 
This applies, e.g., to oxygen-depleted, Pr or Zn-doped 
YBCO (see also Refs. [|jj2|). 

Note that our analysis in this section is based on the 
theory in Ref. |l| (see Eq. @). The core of this theory 
is the introduction of the pair-breaking effect. This effect 
requires singlet pairing and is applicable to both s and 
rf-wave symmetries of the order parameter. There is, of 
course, a quantitative difference between these two cases, 
but the isotope effect for the penetration depth should be 
observed, regardless of the symmetry. In this paper, in 
analogy with Ref. || we make specific calculations for the 
s-wave scenario. Therefore, these calculations are valid, 
even quantitatively for conventional superconductors, as 
well as for some oxides (e.g. Nd-based cuprates, Ba-Pb- 
Bi-O). As for the YBCO compound, our calculations of 
the isotope shift for T c caused also by magnetic impu- 
rities, revealed excellent agreement with the data (see 
Ref. |, Fig. la) which is interesting from the point of 
view of analyzing the complex problem of the symme- 
try. In this paper, we employ a similar approach. It 
is essential, however, to note that, as mentioned above, 
the qualitative picture is similar for both scenarios (s or 
d- wave). 

The results presented in Figs. 1 and 2 were calcu- 
lated on the basis of the parameters for Zn-doped YBCO 
(YBCZnO) and consider solely the influence of magnetic 
impurities on f3. The situation is more complicated if, 
in addition, non-adiabaticity is present in the system. In 
this case f3 m is given by Eqs. ( |l3| ) to (|l7|) and a by 
Eq. (Q). To stress the influence of non-adiabaticity on 
(3 m through ao we write f3 m +na in the following. We have 
calculated the IC of the penetration depth near T c for Pr- 
doped YBCO (YPrBCO) using the parameters derived in 
Ref. | The results are shown in Figs. 1 (dotted and dot- 
dashed lines) and 3 for the temperature and magnetic 
impurity concentration dependence, respectively. 

Several comments can be made by analyzing these fig- 
ures. Fig. 1 shows that there is no qualitative change in 
f3(T) (at fixed T c ) when adding the non-adiabatic contri- 
bution. Non-adiabaticity merely induces a stronger tem- 
perature dependency. This may, however, be specific to 
the materials studied here. On the other hand, contrary 
to the case of magnetic impurities alone (Fig. 2), there 
is a qualitative difference (besides the opposite signs) be- 
tween the IC of T c and S taken as a function of T c and in 
the presence of non-adiabaticity. To understand this, one 
has to come back to Eqs. (p^)-(p^|). Let us first rewrite 
Eq. (|l3|) in the form 

0m = Pm ~ Pa = -~ (fia m - f a ) , (18) 

where fo and /i are defined in Eqs. (p^) , (|l5|) . It turns 
out that /i(T/T c ) and fo{T/T c ) are never more than sev- 



eral percents apart from each other for given T and T c . 
The shape of [3 m+na in the presence of non-adiabaticity 
and magnetic impurities is thus mainly determined by 
the difference a m — ao [where a m is given by Eq. (|lq ) 
and ao = a p h + a na by Eq. (0)]. Using a polynomial 
expansion of T c {n) (as done in Ref. ||) one can show that 
ao saturates with decreasing T c [a na (T c — ► 0) — > 7, with 
7 defined in Eq. (|J)], whereas a m increases as n 2 M in 
the same limit (see Ref. ^|). Thus, at high impurity dop- 
ing, the behaviour of [3 m+na is dominated by the denom- 
inator of Eq. (|l6|), that is by the magnetic impurities. 
In the opposite low doping regime, f3 m + n a is dominated 
by the non-adiabatic part in a m and ao- Accordingly, 
Pm+na behaves similarly near optimal T c and near the 
60K plateau (in YBCO), where dT c /dn is small. Thus, 
as can be seen on Fig. 3b, (3 m +na behaves approximately 
as — a na in these regions but is smaller by an amount ap- 
proximately equivalent to ao. As one gets away from 
the plateaus, the magnetic impurity part of a m — ao 
becomes dominant, leading to a departure from the be- 
haviour of a na . The fact that m +na is mainly deter- 
mined by a m — ao (instead of a m for the IC of T c ) ex- 
plains why the change of curvature of f3 m+na (see Fig. 3b) 
takes place at higher T c 's (lower magnetic impurity con- 
centrations) than the IC of T c , a m (see also Fig. 2a of 
Ref. ||). The change of curvature should be measurable 
since the precision of measurements increases with de- 
creasing dopingB. Accordingly, the effect should also be 
observable for oxygen-depleted YBCO for T c near its op- 
timal value and below but near the 60K plateau. 

We have not calculated the penetration depth IC for 
oxygen-depleted YBCO near T c because the dependen- 
cies F s (ie) (the spin- flip scattering amplitude as a func- 
tion of oxygen doping) and n(x) are not well known for 
this material. This was not a major problem for the 
calculation of the IC a m of T c , because a m is mainly 
determined by the non-adiabatic channel in the regime 
where the experiments were performed and is affected by 
magnetic impurities in the higher doping (i.e. lower T c , 
away from the plateaus) regime (see Eq. (Jl^) and Ref. ||). 
However, it might be of importance for the calculation of 
the IC of 5 because of the more complicated structure of 
the relation p m (T s ). This case is discussed in more detail 
in Ref. I 

Finally, one should note that there is no apparent rea- 
son for the IC of 6 to be restricted to values below 0.5. In 
fact, the calculations show that at high enough dopant 
concentration, the effect should be very large (this is also 
true at T = as shown in the next section) . 



B. Magnetic impurity contribution to (3 at T = 

Since we want to determine solely the contribution of 
non-adiabaticity and magnetic impurities to the IC, one 
is interested in eliminating the presence of the BCS con- 
tribution caused by the factor <p(T/T c ) and the usual 



G 



isotopic dependence of T c . Besides the procedure pre- 
sented in the last section there is another way to avoid 
Pbcs- Indeed, the BCS part becomes negligibly small 
(ip ~ 1 — (T/T c ) 4 ) as one lowers the temperature. At 



T = its contribution is zero. The penetration depth 
calculated by Skalski et al. at zero temperature is given 
by S- 2 = ~(4:Trne 2 /mc 2 )K(uj = 0, q = 0) withN 



£(0,0) 



l + rr 



R(fj) 



Tfj- 



3 V 



7T 

4^ + 1 



(19) 



for r, r\ < 1 (with f(rj) — arcos^) or T < 1, r\ > 1 (with f(rj) = arcoshry) and 



^(0,0) 



i + ry 



l) (T-R(T)) 



R 1 (77) (arcosh?/ — 2artanh7?.) 



(20) 



r-i\ 



for r,?y > 1. We have introduced the notati on T = 

r s /A, 77 = 77F = r 2 /A, 77 = ryr s , i?(x) = ^\i-x 2 \ 

with x = T, 77 and K = [(r-l)(r7- I)/ (V +\){f]+l)] 1/2 . 
A = A(T = 0, T s ) is the order parameter in the presence 
of magnetic impurities. Eqs. (^9|),(|20|) are valid when 
T s <C IV These two scattering amplitudes, T 2 and T s 
(r s = Ti - T 2 ), defined by Abrikosov and Gor'kovta, de- 
scribe the direct and exchange scattering, respectively. 
One can calculate the magnetic impurity contribution to 



the IC at T = from Eq. ([!]) in a straightforward way 
using Eqs. jl9| ) and (20). The result can be written as 



Pm(T = 0) 



a A K\ + K 2 
"~2~ if (0,0) 



(21) 



where a a is defined below, K(0, 0) is given by Eqs. (|lS|), 
© and 



A'i = — 



1 + 3f r 3 (*_m\ ± i 



V 



-3 



#(77) y i?(77) 2 \R(fj) 



V 



•m - 1 



K 2 = ffj- 3 (2 - J77) 

for r,77 < 1 (upper sign, /(77) = arcos77) or T < 1, fj > 1 (lower sign, /(?7) = arcoshfy) and 



#1 = - 



1 + 3r77' 



r - 1 



1 



R(fj) 



may 



i i?(r) i?(77) 

(arcosh77 — 2artanh72.) 



(arcosh?7 — 2artanh72.) 



1 + f 77- 



V 



, f(f-l) 
' i?(f) 3 



7?. 



r. 



r 2 



If' 



K 2 = -[-fj 2 + l)[T- R(T)] - -7]R(T) -n 



K 2 \R(T) 2 
1 



^(4 



R(a) 2 
r-i 



J2(T); 



LiJ(T) 
f 2 /2 



11 "3* 



-—7) 

2 'i?(r) 



a 77 



2?7r 



[f-fl(f)] 
r - 1 



R(T) 



1 



2i?(r) 2 



(22) 



(23) 



for T,t7 > 1 and AT C /T c (r - 1) « 1. The last condi- 
tion expresses the fact that the calculation is not valid 
in the immediate vicinity of f = 1. Eq. (pl| ) contains 
a a which is the IC of the order parameter A. In strong- 

ta (£)-(-Hr+*<r>] + M> 

K 2T 



coupling systems, aA has to be calculated numerically 
using Eliashberg's equations. In the following we calcu- 
late the IC in the framework of the BCS model where c*a 
can be calculated analytically. Indeed, from the relations 

: f < 1 

f , , (24) 
arctani?^)- 1 : T>1 v ' 



derived by Abrikosov and Gor'kov [Aq = A(T = 0,T S = 0) is the order parameter in the absence of magnetic 
impurities] one obtains 
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OiA = ClAo 



(i-if)" 1 



r , R(T) 

1 arctani? r - 1 — 

2 2f 



r < 1 
f > l 



(25) 



In the BCS approximation one further has oa = cxq, 
where the last quantity was defined before as the IC of 
T c q, that is, in the absence of magnetic impurities. As 
for the calculation near T c , the isotope coefficient m at 
T = is independent of any free parameter. Indeed, T s is 
determined by the Abrikosov-Gor'kov equation relating 
T c and T c o, whereas I^ (the direct scattering amplitude) 
is determined via normal state properties of the material. 
Thus the relation j3 m (T c ) is also universal at T = 0. 

Since the penetration depth S does not depend explic- 
itly on T c at T = [see Eqs. @ and @], (3 m {T = 0) is 
independent of the isotope coefficient a. Consequently, 
the contributions of magnetic impurities and the non- 
adiabatic channel to the IC of 8 are simply additive at 
T = 0. This has to be seen in contrast to the calculation 
of the last section performed near T c . In the following we 
neglect the additive non-adiabatic additive contribution 
and show only j3 m . 

The numerical calculation of j3 m from Eqs. (p2|)-(p5|) is 
shown in Figs. 4a and 4b as a function of T c (that is of 
the magnetic impurity concentration) for two values of 
the ratio r) = r 2 /r s (> 1). Fig. 4a displays the case of 
a conventional superconductor and Fig. 4b was obtained 
with parameters corresponding to the situation of high- 
temperature superconductors. The result shows that in 
analogy with the IC of T c in the presence of magnetic 
impurities and with the calculations of the last section, 
\Pm{T = 0)| increases with increasing magnetic impu- 
rity concentration. However, the two coefficients a and 
(3 have opposite signs. 

There are two main differences between the situations 
of conventional and high-temperature superconductors. 
On the one hand, the isotope coefficient j3 m at T = OK of 
high-T c superconductors remains small over a relatively 
large concentration range of magnetic impurities. The ef- 
fect becomes sizeable in the gapless regime. The increase 
of \Pm(T = 0)| is much larger at low doping for conven- 
tional superconductors. As a consequence, on the other 
hand, the effect of the parameter r\ is sizeable only in 
the case of conventional superconductors and in the low 
impurity concentration regime. The calculations suggest 
also that for both types of superconductors the largest ef- 
fect is observed at low T c where the superconductor is in 
the gapless state (T, fj > 1). The values of the IC are then 
of the same order as those of the IC of T c . It would be 
interesting to perform experiments both on conventional 
and high-temperature superconductors. As discussed be- 
fore, Zn- or Pr-doped YBCO are interesting candidates 
since the IC of T c as a function of doping (which is nec- 
essary in order to csdailate /3 tot = f3 ph + (3 m + na ) is 
known in both cases BoH 



IV. ISOTOPE EFFECT OF S IN PROXIMITY 
SYSTEMS 

In this last section we consider another physical situa- 
tion, a proximity system, in which a factor not related to 
lattice dynamics, induces an IE of the penetration depth. 
We have shown in Ref. |^ that the presence of a normal 
film on a superconductor induces a change in the IC of 
T c . Here we show that it also induces an isotopic shift of 
(5. We study a proximity system of the type S — N, where 
S is a weak-coupling superconductor and N is a metal or 
a semiconductor (e.g. Pb-Ag). Let us note T c0 the value 
of the critical temperature for an isolated S film. As is 
known, the proximity effect leads to a critical tempera- 
ture T c that differs from T c0 (see, e.g., Refs. p|,p3|). The 
proximity effect also affects the shielding of a weak mag- 
netic field. The most dramatic effect of the normal layer 
on the penetration depth is seen in the low temperature 
regime (T/T c < 0.3). Although the penetration depth of 
a pure conventional superconductor is only weakly tem- 
perature dependent in this regime (8 ~ \J\ — T A /T^ ) 
the presence of the normal layer induces a temperature 
dependence through the proximity effect. 

Consider a proximity system S — N and assume that 
8 < Lpj *C £jv (which is certainly satisfied in the low- 
temperature regime since £jv = Tlvf/2-kT; see, e.g., 
Ref. HI) , where L jy and £jv are the thickness and the co- 
herence length of the normal film, respectively (cf. Ref. ^ 
and references therein). In this case, the peneipation 
depth 8 has been calculated by one of the authorsc-3: 

8- 3 =a N <f> , (26) 

where a at is a constant depending only on the material 
properties of the normal film (it is independent of the 
ionic mass) and 

$ = ^E 2 2,1 ' Pn = l+6t^lTT . (27) 

x n = u n /es(T) with w„ = (2n + l)7rT (the Mastubara 
frequencies) and es(T) is the superconducting energy 
gap of S. In the weak-coupling limit considered here, 
es(0) = £ttT c q with e ~ 0.56. The dimensionless pa- 
rameter t — £/Sq with £ = Ln/Lq and-.S'o = ro/7rT c o 
(To ~ 1/Lq is the McMillan parametertj) . Lat and Lq 
are the thickness of the normal film and some arbitrary 
thickness, respectively (in the following we take Lq = Lg, 
the thickness of the superconducting film) . 

Since 8 depends non-linearly on T c q (through the prox- 
imity parameter t), the penetration depth will display an 
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isotope shift due to the proximity effect. From Eq. ( p6| ) 
one can calculate the isotope coefficient f3 pr0 x [as defined 
by Eq. (|])] of the penetration depth due to the proximity 
effect: 

3$ „>o <pI + 1 V pW x * + 1 J 

where ao is the IC of T c q for the superconducting film S 
alone. 

This result has several interesting properties. First, 
it turns out that the isotope coefficients f3 prox and ao 
have generally opposite signs (see below). Secondly, the 
IC depends on the proximity parameter £, that is on the 
thickness ratio I — Ljq jLg of the normal and supercon- 
ducting films and on the MacMillan tunneling parameter 

r . 

The dependency on the ratio £ has. already been men- 
tioned in connection with the IE of Tj3'U. Fig. 5 shows the 
isotope coefficient (3 prox as a function of the parameter t 
for fixed temperatures. \l3 prox \ decreases with increasing 
t. There are two ways to change the value of t ~ £/Sq. 
One of them is an increase of the parameter £ = Ljq/ L$ 
(e.g. an increase of the thickness of the normal film Ljy). 
Note that there is a lower bound £ m i n to the value of_£, 
since it was assumed that 5 < Ln to obtain Eq. (^)lj. 
The IC \(3 prox \ decreases with increasing thickness ra- 
tio (Fig. 5 then corresponds to the values So = 0.2 for 
£ G]^mm,l])- However, the higher the temperature the 
smaller the decrease of (3 pr0 x which then rapidly reaches 
saturation with increasing thickness. 

Instead of changing the thickness ratio £ one can vary 
the parameter Tq. One possibility is to use different films, 
another is to modify the quality of the S — N interface. 
Fig. 5 then implies that the \/3 pr0 x\ decreases with de- 
creasing tunneling parameter To- It would be interesting 
to perform an experiment where the isotope effect of T c 
and 6 are determined for a proximity system with differ- 
ent values of the thickness ratio £. 

Note that the result (^8|) has another interesting fea- 
ture ; one can see that the isotopic shift of the penetra- 
tion depth is temperature dependent. Fig. 6 shows the 
temperature dependence of (i prox for two sets of param- 
eters. The two upper curves are obtained for So = 0.2 
and two thickness ratios £ = 0.5, 1; the two lower curves 
are for So = 1 and £ — 0.5, 1. From Fig. 6 one concludes 
that |/3 P roa;| increases with increasing temperature. Fur- 
thermore, a weak tunneling between the superconducting 
and the normal film (smaller So) induces a strong temper- 
ature dependence. Finally, the temperature range over 
which the IC remains unchanged is broader for a larger 
than a smaller tunneling parameter. 

Therefore, the presence of a normal layer on a su- 
perconductor induces an isotope shift of the penetration 
depth. Moreover, in a proximity system the IC depends 
on the thickness ratio and the temperature. The values of 
the isotope coefficient are such that they are measurable 
in conventional and high-temperature superconductors. 



V. CONCLUSION 



We have shown that the value of such a fundamental 
parameter as the penetration depth S depends on iso- 
topic substitution. This unconventional dependency is 
caused by number of factors not related to the pairing 
mechanism. 

At first, we focused on the non-adiabatic IE. This ef- 
fect, manifested in the dependence of the charge-carrier 
concentration, on the ionic masses for doped systems such 
as cupratesou, leads to a noticeable isotopic shift of the 
penetration depth 5. For this case, we established a gen- 
eral relation between the isotope coefficient of T c and of 
S for London superconductors (see Eq. (|l0|)). The iso- 
tope coefficient of <5 is a constant independent of doping 
or temperature; we determined its value for oxygen de- 
pleted and Pr substituted YBa2Cu307_<5 as well as for 
La2- 2: Sr 2 ,CuO | +_ Note that similar effects are expected 
in manganitesEHl where the charge-transfer processes also 
involve the non-adiabatic oxygen ions. 

The presence of magnetic impurities also leads to a 
non-trivial isotope effect for S. We studied this effect for 
temperatures near T c and at T = 0. It turns out that the 
coefficient obeys a universal law; it is determined by the 
experimental values of T c in the absence and presence of 
magnetic impurities [Eqs. ( |l3| ) and fl2l])] as well as the 
isotope coefficient of T c in the absence of magnetic im- 
purities. The combined effect of magnetic impurities and 
non-adiabaticity was also discussed. Finally, we showed 
that the penetration depth for proximity systems also 
display an isotopic shift. It is interesting that in the case 
of a proximity system or in the presence of magnetic im- 
purities the IC for d is temperature dependent. 

All these effects should be measurable in conventional 
and in high-temperature superconductors and several ex- 
periments were proposed. In the case of high-T c mate- 
rials it would be of interest to study oxygen depleted 
YBa2Cu307_a as well as the Pr- and Zn-substituted 
systems, since the IE of T c has been successfully de- 
scribed in these materials (see Refs. |l|-||). The case of 
Zn-substituted YBa2Cu30g+x is here of special inter- 
est, since experiments show that the depression of T c is 
mainly due to the change in magnetic impurity concen- 
tration and our theory is free of any adjustable parameter 
in this case. 

Note, finally, that not only T c and 5 exhibit the uncon- 
ventional isotope effects studied here and in Refs. 
but any quantity that depends on the charge carrier 
density of the superconducting condensate or/and on T c 
(e.g. the specific heat). 
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Figure captions: 

Fig. 1: Temperature dependency (near T c ) of the isotope coeff. j3 a) in the presence of magnetic impurities only, 
with T c /T c0 ~ 0.35 (solid), 0.5 (dashed line), and a = 0.025 (corresponds to YBa2(Cui_ x Zn x )30 7 _,5); b) 
in the presence of magnetic impurities and non-adiabaticity for T c /T c q ~ 0.5 (dash-dotted), 0.6 (dotted line). 
«o = 0.025, 7 S = 0.16 (corresponds to Yi_ a; Pr :I ;Ba2Cu307_i, see text). 

Fig. 2: Isotope coefficient (3 m near T c as a function of magnetic impurity concentration for T/T c = 0.75 (solid line), 
0.85 (dashed) and 0.95 (dotted). a = 0.025 (corresponds to the situation of YBa^Cui-zZn^C^-^). 

Fig. 3a: (3 m +na as a function of magnetic impurity concentration in the presence of non-adiabaticity for T/T c = 
0.75 (solid line), 0.85 (dashed) 0.95 (dotted). a = 0.025, 7 = 0.16 and f = 123K (parameters for 
Y 1 _,Pr ;c Ba 2 Cu30 7 _5). 

Fig. 3b: Low doping regime of Fig. 3a. In addition, the case without non-adiabaticity is added for comparison (long 
dashed line, corresponding to YBa^Cui-^Zn^C^-a with T/T c = 0.9). 

Fig. 4a: Isotope coeff. (3 m for conventional superconductors at T = as a function of magnetic impurity concentration 
for T 2 /T s = 10 (solid line) and 50 (dashed). a = 0.3. 

Fig. 4b: Isotope coeff. f3 m for high-T c superconductors at T — as a function of magnetic impurity concentration 
for r 2 /r s = 10 (solid line) and 50 (dashed) a = 0.025. 

Fig. 5: Isotope coeff. f} prox for a proximity system as a function of t for T/T c = 0.1 (solid line), 0.2 (dashed) and 0.3 
(dotted). a a = 0.5. 

Fig. 6: Isotope coeff. (3 prox for a proximity system as a function of T/T c for So = 0.2: I = 1 (dashed) I = 0.5 (solid) 
and So = 1: 1 = 1 (small dots) I — 0.5 (dotted). a = 0.5. 
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